正则变换 生成函数

1. 正则变换

正则变换:保证哈密顿方程变换前后方程形式不变. 引入泊松括号表示: \[[\tilde{f}, \tilde{g}]_{QP} = [f, g]_{qp}.\]

1.1. 正则变换的定义-–—泊松括号条件

传统相空间中 \(q, p = q_1, q_2, \dots, q_D, p_1, \dots, p_D\). 变换方程:

\begin{cases} Q_k = Q_k(q, p) = Q_k(q_1,\dots, q_D, p1, \dots, p_D)\\ P_k = P_k(q, p) = P_k(q_1,\dots, q_D, p1, \dots, p_D) \end{cases}

构成正则变换的充要条件为: \[[Q_k, Q_l]_{qp} = 0,\;[Q_k, P_l]_{qp} = \delta_{kl},\;[P_k, P_l]_{qp} = 0.\] 由于

\begin{aligned} [Q_k, P_l]_{qp} &= \left[ \frac{\partial Q_k(q, p)}{\partial q_j}\frac{\partial P_l(q, p)}{\partial p^j}- \frac{\partial Q_k(q, p)}{\partial p_j}\frac{\partial P_l(q, p)}{\partial q^j}\right] \\ &=a_{kj}b_l^j = \left( ab^T \right)_{kl}. \end{aligned}

故得: \(\left( ab^T \right)_{kl} = \delta_{kl}\), 即 \(ab^T = I_D\) 为单位阵时, 变换为正则变换.

1.2. 辛坐标(symplectic coordinates) 辛矩阵

1.2.1. 相空间辛坐标

\[\gamma_k = q_k, \gamma_{D+k} = p_k,\;k = 1, 2, \dots, D\]. 变换 \(Q_k(q, p), P_k(q, p)\) 转化为 \(\Gamma_i = \Gamma_i(\gamma)\).

雅可比矩阵
\[J_{ij} = \frac{\partial \Gamma_i(\gamma)}{\partial\gamma_j}, \;J_{2D\times 2D} = \begin{bmatrix}\partial_q Q & \partial_p Q \\ \partial_q P & \partial_p P\end{bmatrix}\]

有 \[[\gamma] = \begin{pmatrix}\gamma_1\\ \vdots\\ \gamma_{2D}\end{pmatrix}, \;[\dot{\gamma}] = \begin{pmatrix}[\dot{q}]\\ [\dot{p}]\end{pmatrix}.\] 相空间函数写成: \(f(q, p) = f(\gamma)\), 有: \[\left[ \frac{\partial f}{\partial \gamma} \right] = \begin{pmatrix} \left[ \dfrac{\partial f}{\partial q} \right] \\ \left[ \dfrac{\partial f}{\partial p} \right]\end{pmatrix}.\] 同样有: \[ [\dot{\Gamma}] = \begin{pmatrix} [\dot{Q}] \\ [\dot{P}] \end{pmatrix}, \; \left[ \frac{\partial f}{\partial \Gamma} \right] = \begin{pmatrix} \frac{\partial f(\Gamma)}{\partial \Gamma_1}\\ \vdots \\ \frac{\partial f(\Gamma)}{\partial \Gamma_{2D}}\end{pmatrix}.\]

因为 \([\dot{\Gamma}] = J \dot{\gamma}\), 有: \[ \left[ \frac{\partial f}{\partial \gamma} \right] = J^T \left[ \frac{\partial f}{\partial \Gamma} \right]. \]

1.2.2. 辛矩阵

\[ s_{2D\times 2D} = \begin{bmatrix} 0 & I_D \\ -I_D & 0\end{bmatrix}, \;s_{ij} = \delta_{i,j-D} - \delta_{i,j+D}.\] 易证:

  • \(s^T = -s\),
  • \(ss = -I_{2D}\),
  • \(s^T s = I_{2D}\),
  • \(\det(s) = 1\).

1.3. 辛形式下的方程

  • 哈密顿方程: \[ [\dot{\gamma}] = s \left[ \frac{\partial H}{\partial \gamma} \right]. \]
  • 泊松括号: \[[f, g]_{qp} = \left[ \frac{\partial f(\gamma)}{\partial \gamma} \right]^T s \left[ \frac{\partial g(\gamma)}{\partial\gamma} \right].\]

1.4. 正则变换的泊松括号条件的辛形式: \(JsJ^T = s\)

由前述, 正则变换满足: \[[Q_k, Q_l]_{qp} = 0,\;[Q_k, P_l]_{qp} = \delta_{kl},\;[P_k, P_l]_{qp} = 0.\] 即 \[ [\Gamma_k, \Gamma_l]_{qp} = s_{kl}. \] 辛形式下的泊松括号有: \[ [\Gamma_k, \Gamma_l]_{qp} = \left[ \frac{\partial \Gamma_k}{\partial \gamma} \right]^T s \left[ \frac{\partial \gamma_l}{\partial \gamma} \right] = \left[ \frac{\partial \Gamma_k}{\partial \gamma_i} \right]^T s_j^i \left[ \frac{\partial \gamma_l}{\partial \gamma^j} \right] = \left( JsJ^T \right)^T_{kl}\] 故得: \[JsJ^T = s.\]

1.5. 正则变换的可逆性

由前述, 有: \(\left| JsJ^T \right| = |s|\), 可得 \({|J|}^2 = 1 \ne 0\). 故\(J\)可逆, 易得逆矩阵 \(J^{-1}\) 有: \[J^{-1} = \begin{bmatrix} \partial_Q q & \partial_P q \\ \partial_Q p & \partial_P p \end{bmatrix}\]

1.6. 正则变换的直接条件: \(J^{-1} = -sJ^T s\)

存在 \(J^{-1}\) 使得 \(J^{-1} = -sJ^Ts\), 则变换为正则变换. 此时有: \[\left(\frac{\partial q}{\partial Q}\right)=\left(\frac{\partial P}{\partial p}\right)^T, \;\left(\frac{\partial q}{\partial P}\right)=-\left(\frac{\partial Q}{\partial p}\right)^T, \;\left(\frac{\partial p}{\partial Q}\right)=-\left(\frac{\partial P}{\partial q}\right)^T, \;\left(\frac{\partial p}{\partial P}\right)=\left(\frac{\partial Q}{\partial q}\right)^T,\]

1.7. 正则变换的拉格朗日括号判据: \(J^T s J = s\)

从泊松括号出发: \(JsJ^T = s\), 且 \(ss = -I_n\).

  1. 左乘 \(sJ^{-1}\): \((sJ^{-1})JsJ^T = (sJ^{-1})s\Rightarrow -J^T = sJ^{-1}s\)
  2. 右乘 \(sJ\): \(J^Tsj = -sJ^{-1}ssJ = s\).

Vice versa. 充要条件.

由此得到正则变换三大判据:

  1. 泊松括号判据: \(JsJ^T = s\),
  2. 直接判据: \(J^{-1} = -sJ^Ts\),
  3. 拉格朗日括号判据: \(J^TsJ = s\).

1.7.1. 拉格朗日括号的长形式

\[\left\{ \gamma_i, \gamma_j \right\}_{QP} = \left[ \frac{\partial Q_k}{\partial \gamma_i}\frac{\partial P^k}{\partial\gamma_j}- \frac{\partial Q_k}{\partial\gamma_j}\frac{\partial P^k}{\partial\gamma_i}\right].\] 正则变换时, 有: \(\{\gamma_i, \gamma_j\}_{QP} = s_{ki}\).

泊松括号与拉格朗日括号的关系:

  • 泊松括号为确定分子指标,对分母指标求和.
  • 拉格朗日括号为确定分母指标,对分子求和.

1.8. 正规群

群:定义了乘法且满足以下性质的集合.

  1. 封闭性:两个元素之积仍在集合内.
  2. 存在单位元:对应于作用前后不做任何改变.
  3. 存在逆元素:一个元素与其逆元素之积为单位元素.
  4. 乘法满足结合律(但不一定满足交换律).

正则变换构成一个群

  1. 封闭性: 若 \(A, B\) 为正则变换, 则 \(C = AB\) (先做 \(B\) 变换, 再做\(A\)变换) 也是正则变换.
  2. 存在单位元 \(Q_k = q_k,\;P_k = p_k,\;J = I_n\), 显然 \(I_n sI_n = s\).
  3. 存在逆元素:正则变换可逆,逆变换也为正则变换.
  4. 乘法结合律:由矩阵乘法易证.

1.9. 泊松括号形式不变性

泊松括号具有形式不变性 \([f, g]_{QP} = [f, q]_{qp}\), 等价于相空间变换为正则变换 \(\left[ \frac{\partial f}{\partial\gamma} \right] = J^T \left[ \frac{\partial f}{\partial\Gamma} \right]\) . 以下给出证明:

  • 形式不变性⇒ 正则变换: 取 \(f = Q_i\), \(g = Q_j\), 易得 \[[Q_i, Q_j]_{QP} = [Q_i, Q_j]_{qp} = 0.\] 同理可得构成正则变换的其余充要条件, 故成立.
  • 正则变换⇒ 形式不变性: \[[f, g]_{qp} = \left[ \frac{\partial f}{\partial \gamma} \right]^T s \left[ \frac{\partial g}{\partial\gamma} \right], \;\left[ \frac{\partial f}{\partial\gamma} \right] = J^T \left[ \frac{\partial f}{\partial\Gamma} \right].\] 可得: \[[f, g]_{qp} = \left[ \frac{\partial f}{\partial \Gamma} \right]^T JsJ^T \left[ \frac{\partial g}{\partial\Gamma} \right],\;JsJ^T = s.\] 即: \[[f, g]_{qp} = \left[ \frac{\partial f}{\partial \Gamma} \right]^T s \left[ \frac{\partial g}{\partial\Gamma} \right] = [f, g]_{QP}.\blacksquare\]

1.10. 哈密顿方程的形式不变性

在 \(q, p\)-系中, 成立: \[\dot{q}_k = \frac{\partial h(q, p, t)}{\partial p_k}, \;\dot{p}_k = -\frac{\partial h(q, p, t)}{\partial q_k}.\] 即 \[[\dot{\gamma}] = s \left[ \frac{\partial h}{\partial \gamma} \right].\] 有正则变换

\begin{cases} q_k = q_k(Q, P)\\ p_k = p_k(Q, P) \end{cases}

则哈密顿量有 \(H(Q, P, t) = h \left[ q(Q, P), p(Q, P), t \right]\). 在 \(Q, P\)-系中, 哈密顿方程有相同形式: \[\dot{Q}_k = \frac{\partial H(Q, P, t)}{\partial P_k}, \;\dot{P}_k = -\frac{\partial H(Q, P, t)}{\partial P_k}.\]

以下给出证明. 利用变换辛形式: \[\left[\dot{\Gamma}\right] = J[\dot{\gamma}], \;\left[ \frac{\partial f}{\partial\gamma} \right] = J^T \left[ \frac{\partial f}{\partial\Gamma} \right].\] 且 \([\dot{\gamma}] = s \left[ \partial_{\gamma}h \right]\), 可得: \[\left[ \dot{\Gamma} \right] = JsJ^T \left[ \frac{\partial H}{\partial\gamma} \right] = s \left[ \frac{\partial H}{\partial\gamma} \right].\blacksquare\]

注意当正则变换显含时间时, \(\left[ \dot{\Gamma} \right]=J[\dot{\gamma}]\) 不成立. 此时有: \[\frac{\mathrm{d}\Gamma_k}{\mathrm{d}t} = \frac{\partial\Gamma_k}{\partial t} + \frac{\partial\Gamma_k}{\partial\gamma_l}\dot{\gamma}^l.\] 当变换显含时间且时间不参与变换时,哈密顿量需要修正. \[\dot{Q}_k(q, p, t) = [Q_k, h]_{qp} + \frac{\partial Q_k(q, p, t)}{\partial t}.\] 设 \(H(q, p, t) = h(q, p, t) + g(q, p, t)\), \(H(Q, P, t) = H[q(Q, P), p(Q,P), t]\), 此时有:

\begin{aligned} \dot{Q}_k &= [Q_k, h]_{qp} + \frac{\partial Q_k}{\partial t} \\ &= \left[ Q_k, h \right]_{qp} + [Q_k, g]_{qp}\\ &= [Q_k, H]_{qp}=[Q_k, H]_{QP} = \frac{\partial H(Q, P, t)}{\partial P_k}. \end{aligned}

同理可得: \[ \frac{\partial P_k(q, p, t)}{\partial t} = [P_k, g]_{qp}, \;\dot{P}_k = -\frac{\partial H(Q, P, t)}{\partial Q_k}.\] 故: 变换显含时间时, 需要 \(H = h + g\), 其中 \[ g = \begin{cases}\frac{\partial Q_k}{\partial t} = [Q_k, g]_{qp}\\ \frac{\partial P_k}{\partial t} = [P_k, g]_{qp}\end{cases}. \]

扩展理论:扩展理论将时间\(t\)作为变换坐标与其他坐标一起变换.

2. 生成函数

2.1. 正则变换的又一种定义

由变分原理知,使得作用量积分取驻值的相空间路线满足哈密顿方程. \[\int_{t(1)}^{t(2)}\left[ p_k\dot{q}^k - H(q, p, t)\right]\mathrm{d}t\stackrel{Q,\;P}{\Longrightarrow} \int_{t(1)}^{t(2)}\left[ P_k\dot{Q}^k - K(Q, P, t) \right]\mathrm{d}t.\]

正则变换(变分原理判据;条件;定义)
满足下式的可逆变换 \(Q_k = Q_k(q, p, t)\), \(P_k = P_k(q, p, t)\) 称为正则变换: \[P_k\dot{Q}^k - K(Q, P, t) = p_k\dot{q}^k - H(q, p, t) - \frac{\mathrm{d}U}{\mathrm{d}t}.\] 即 \[p_k\mathrm{d}q^k - P_k\mathrm{d}Q^k + (K - H)\mathrm{d}t = \mathrm{d}U.\] 即上式左式是某个多元函数 \(U\) 的全微分. 此定义与泊松括号定义一致. \(U\) 的含义即下一节所述生成函数.

2.2. 生成函数

由前述, 函数 \(U(q, Q, t)\) 决定了正则变换, 称为正则变换的生成函数. 其偏微分有: \[p_k = \frac{\partial U}{\partial q_k},\;P_k =-\frac{\partial U}{\partial Q_k}, \;K - H = \frac{\partial U}{\partial t}.\] 变换可逆, 可求得 \(q_k, p_k\), \(Q_k, P_k\) 关系. \[K(Q, P, t) = H[q(Q, P), p(Q, P), t] + \frac{\partial U(q, Q, t)}{\partial t}\Big{|}_{q_k = q_k(Q, P, t)}.\] 只有\(U\)不显含时间时,哈密顿量不变.

第一类生成函数
记上述 \(U(q, Q, t)\) 为 \(U_1(q, Q, t)\) 即第一类生成函数.
第二类生成函数
\(U_2(q, P, t) = U_1 + P_kQ^k\). \[ p_k = \frac{\partial U_2}{\partial q_k},\;Q_k = \frac{\partial U_2}{\partial P_k}, \;K = \frac{\partial U_2}{\partial t} + H. \]
第三类生成函数
\(U_3(p, Q, t) = U_1(q, Q, t) - p_kq^k\). \[q_k = -\frac{\partial U_3}{\partial p_k},\;P_k = -\frac{\partial U_3}{\partial Q_k}, \;K = \frac{\partial U_3}{\partial t} + H.\]
第四类生成函数
\(U_4(p, P, t) = U_3+P_kQ^k\). \[q_k = -\frac{\partial U_4}{\partial p_k},\;Q_k = \frac{\partial U_4}{\partial P_k}, \;K = \frac{\partial U_4}{\partial t} + H.\]

除了这四种生成函数,还有混合型生成函数.

2.2.1. 记忆方法

four-gen-func.svg
图1  \(U_i\) 所在象限决定了 \(p\mathrm{d}q\) 和 \(P\mathrm{d}Q\) 的符号.
  • \(\mathrm{d}U_1 = p_k\mathrm{d}q^k - P_k\mathrm{d}Q^k + (K-H)\mathrm{d}t\).
  • \(\mathrm{d}U_2 = p_k\mathrm{d}q^k + Q_k\mathrm{d}P^k + (K-H)\mathrm{d}t\).
  • \(\mathrm{d}U_3 = - q_k\mathrm{d}p^k - P_k\mathrm{d}Q^k + (K-H)\mathrm{d}t\).
  • \(\mathrm{d}U_4 = - q_k\mathrm{d}p^k + Q_k\mathrm{d}P^k + (K-H)\mathrm{d}t\).

2.3. 无穷小正则变换(微分变换)

2.3.1. 无穷小生成元

设有第二类生成函数: \[U_2(q, P, t) = q_kP^k\] 其对应于恒等变换

\begin{cases} Q_k = \partial_{P_k}U_2 = q_k,\\ p_k = \partial_{q_k}U_2 = P_k. \end{cases}

无穷小变换基于恒等变换,其生成函数为: \[U(q, P, t) = q_kP^k + \epsilon G(q, P, t),\;\epsilon\rightarrow 0.\] 有:

\begin{cases} p_k = \frac{\partial U}{\partial q_k} = P_k + \epsilon \frac{\partial G(q, P, t)}{\partial q_k} = P_k + \epsilon \frac{\partial G(q, p, t)}{\partial q_k} + o(\epsilon),\\ Q_k = \frac{\partial U}{\partial P_k} = q_k + \epsilon\frac{\partial G(q, P, t)}{\partial P_k} = q_k + \epsilon \frac{\partial G(q, p, t)}{\partial p_k} + o(\epsilon). \end{cases}

略去高阶小量 \(o(\epsilon)\), 有:

\begin{cases} Q_k - q_k = \epsilon \dfrac{\partial G(q, p, t)}{\partial p_k}, \\ P_k - p_k = - \epsilon \dfrac{\partial G(q, p, t)}{\partial q_k}, \\ K - H = \epsilon \dfrac{\partial G(q, p, t)}{\partial t}. \end{cases}

泊松括号形式为: \[\delta q_k = Q_k - q_k = \epsilon[q_k, G]_{qp},\;\partial p_k = P_k -p_k = \epsilon[G, p]_{qp}.\]

2.3.2. 主动变换

将正则变换看作相点变动—相坐标系不变,相点变动导致相坐标改变.

2.3.3. 运动方程的无穷小正则变换表示

设哈密顿量 \(H(q, p, t) = G(q, p, t)\), 由哈密顿方程可得: \[\delta q_k = \dot{q}_k\mathrm{d}t,\;\delta p_k = \dot{p}_k\mathrm{d}t, \;K = H + \frac{\partial H(q, p, t)}{\partial t}\mathrm{d}t.\] 正则变换有:

\begin{cases} Q_k(t) = q_k(t) + \partial q_k = q_k(t) + \dot{q}_k \mathrm{d}t = q_k(t + \mathrm{d}t), \\ P_k(t) = p_k(t) + \partial p_k = p_k(t) + \dot{p}_k \mathrm{d}t = p_k(t + \mathrm{d}t). \end{cases}

能量有: \[K(Q,P,t) = H(q,p,t) + \dot{H}\mathrm{d}t = H(q,p,t +\mathrm{d}t).\]

  • 等式左侧为经过无穷小变换得到的新相空间坐标.
  • 右侧为运动方程决定的相点在同一相空间新坐标.

⇒ 力学系统在\(t\)到\(t+\mathrm{d}t\)的演化是一个哈密顿量为无穷小生成元的无穷小正则变换.

即力学系统在有限时间间隔的演化等价于一个由以哈密顿为生成元的一系列无穷小正则变换之积构成的正则变换.

2.3.4. 诺特定理的相空间版

一般物理量 \(u(q, p, t)\) 在生成元 \(G(q, p, t)\) 生成的正则变换下变为:

\begin{aligned} \label{orgd9e3e49} \partial u &\equiv u(Q, P, t) - u(q, p, t)\\ &= \left( \frac{\partial u}{\partial q_k}\delta q^k + \frac{\partial u}{\partial p_k}\delta p_k \right) = \epsilon [u, G]_{qp}. \end{aligned}

即 \(u(Q, P, t) - u(q, p, t) = \epsilon [u, G]_{qp}\).

设有不显含时间的物理量 \(G(q, p)\). \[\frac{\mathrm{d}G}{\mathrm{d}t} = [G, H]_{qp} + \frac{\partial G}{\partial t} = [G, H]_{qp}. \] 以\(G\)为生成元,构造无穷小正则变换,由上式知: \[\delta H = H(Q, P, t) - H(q, p, t) = \epsilon [H, G]_{qp}.\] 若这个变换下, 哈密顿量不变, 即 \(\delta H = 0\), 则有 \([H, G]_{qp} = 0\). 可得: \[\frac{\mathrm{d}G}{\mathrm{d}t} = [G, H]_{qp} = -[H, G]_{qp} = 0 \Rightarrow \frac{\mathrm{d}G}{\mathrm{d}t} = 0.\] 作为物理量的 \(G(q, p)\) 是 守恒量 (运动常数). 这就是相空间上的 Noether 定理. 注意 \(G(q, p)\) 不显含时间.

2.3.5. 有限正则变换与 Lie 级数

2.3.6. Liouville 定理

相空间体积 \[\mathcal{V}_t = \int \mathrm{d}q_1(t)\mathrm{d}q_2(t)\cdots \mathrm{d}q_D(t) \mathrm{d}p_1(t)\mathrm{d}p_2(t)\cdots \mathrm{d}p_D(t).\]

Liouville 定理
哈密顿系统的时间演化保持相空间体积不变. 即 \(\mathcal{V}_t = \mathcal{V}_{t_0}\).

已知力学在有限时间间隔的演化等价于以哈密顿量为生成元的一系列无穷小正则变换之积构成的正则变换.

3. 哈-雅理论

通过正则变换求解哈密顿方程.

3.1. 哈-雅方程

无论哪一种生成函数,哈密顿量总是以下列形式变化: \[K = H + \frac{\partial U}{\partial t}.\] 做正则变换使得: \[K = H + \frac{\partial U}{\partial t} = 0.\] 则有:

\begin{aligned} Q_k &= \frac{\partial K}{\partial P_k} = 0 \Rightarrow Q_k = \beta_k = \text{const},\\ \dot{P}_k &=-\frac{\partial K}{\partial P_k} = 0 \Rightarrow P_k = \alpha_k = \text{const}. \end{aligned}

取第二类生成函数 \(U(q, P, t)\). \[p_k = \frac{\partial U_2}{\partial q_k},\;Q_k = \frac{\partial U_2}{\partial P_k}.\] 可得:

\[K = H(q_1, q_2,\dots, q_D, \frac{\partial U_2}{\partial q_1}, \frac{\partial U_2}{\partial q_2},\dots \frac{\partial U_2}{\partial q_D}, t) + \frac{\partial U_2}{\partial t} = 0.\] 上式称为 哈密顿-雅克比方程, 可用于求解生成函数 \(U_2(q, P, t)\). 哈-雅方程的解称为哈密顿作用量函数. (也称哈密顿主函数)

  • 哈-雅方程为一阶非线性偏微分方程, 有 \(D+1\) 个自变量 \(q_1, q_2, \dots, q_D, t\).
  • 哈-雅方程的解含有 \(D+1\) 个积分常数, 最后一个常数不重要, 直接略去.

故方程解含有 \(D\) 个积分常数: \(\alpha_k,\;k=1, 2, 3\dots,D\). 认定 \(P_k = \alpha_k\) (注意到 \(P_k\) 为常数), 即可构造生成函数 \(U_2(q, P, t)\).

归纳出哈雅方程求解方法:

  1. \(p_j = \frac{\partial U_2}{\partial q_j}\) 代入哈密顿量, 得到哈-雅方程: \[K = H \left( q, \frac{\partial U_2}{\partial q}, t \right) + \frac{\partial U_2}{\partial t} = 0.\]
  2. 求解 \(U_2(q1,q2\dots,q_D, \alpha_1,\dots \alpha_D, t)\), 认定积分常数 \(\alpha_j = P_j,\;j=1, 2, \dots, D\). 偏微分方程仅对简单问题可解. (3.2法)
  3. 利用生成函数求 \(Q\): \[Q_j = \frac{\partial U_2(q, \alpha, t)}{\partial \alpha_j} = \beta_j = \text{const}.\] 进而计算出 \(q_j(\alpha, \beta, t),\;j=1, 2, \dots, D\).
  4. 计算 \(p_j = \frac{\partial U_2(q, \alpha, t)}{\partial q_j}\).
  5. 由初始条件确定 \(\alpha, \beta\).

3.1.1. 哈-雅方程的解哈密顿主函数 \(S(q, P, t)\) 的物理意义.

\[\frac{\mathrm{d}S(q, P, t)}{\mathrm{d}t} = \left( \frac{\partial S}{\partial q_k}\dot{q}^k + \frac{\partial S}{\partial P_k}\dot{P}^k \right) + \frac{\partial S}{\partial t}.\] 且 \(p_k = \frac{\partial S}{\partial q_k},\;Q_k = \frac{\partial S}{\partial P_k}, \;\frac{\partial S}{\partial t} = K -H\), \(K = 0,\;\dot{P}_k = 0,\;\dot{Q}_k = 0\). 可得 \[\frac{\mathrm{d}S}{\mathrm{d}t} = p_k\dot{q}^k - H = L.\] 即拉格朗日量. 对上式积分得到: \[S = \int L\mathrm{d}t + \text{const}.\]

  • 作用量的定积分可用于导出哈密顿运动方程.
  • 作用量的不定积分可用于生成正则变换,求解运动方程.

3.2. 分离变量

求偏微分方程特解最常用的方法: 分离变量. 这里将变量分离成和的形式,把循环坐标(可遗坐标)变量分离出来.

3.2.1. 可遗坐标 (循环坐标) ignorable,cyclic coordinate

设 \(q_D\) 为可遗坐标, \(\dot{p}_D = \frac{\partial H}{\partial q_D} = 0\), 即哈密顿量不显含 \(q_D\).

哈-雅方程出发: \[K = H \left( q, \frac{\partial U_2}{\partial q}, t \right) + \frac{\partial U_2}{\delta t} = 0.\] 将哈密顿主函数中 \(S(q, P, t)\) 与可遗坐标 \(q_D\) 相应的函数分离出来 \(S = \alpha_D q_D + \tilde{S}(q_1, \dots, q_{D-1}, t)\). (考虑到 \(P\) 为常数) 代入哈-雅方程, 得: \[H \left( q_1, \dots, q_{D-1},\frac{\partial \tilde{S}}{\partial q_1},\dots, \frac{\partial \tilde{S}}{\partial q_{D-1}},\alpha_D, t \right) + \frac{\partial\tilde{S}}{\partial t} = 0.\] 可得 \[p_D = \frac{\partial S}{\partial q_D} = \alpha_D.\] 故涉及可遗坐标的变换可以是恒等变换: \(U_2 = q_DP_D = q_D\alpha_D\). 存在其他循环坐标同理,可求得特解.

3.2.2. 时间变量

当哈密顿量不显含时间时,哈-雅方程退化为: \[H \left( q, \frac{\partial U}{\partial q} \right) + \frac{\partial U_2}{\partial t} = 0.\] 分离时间变量: 哈密顿主函数 \(S = W(q_1, \dots, q_D) - E t\). 称 \(W(q)\) 为 哈密顿特征函数, 上式代入哈-雅方程: \[H \left( q_1, q_2, \dots, q_D, \frac{\partial W}{\partial q_1}, \frac{\partial W}{\partial q_2},\dots, \frac{\partial W}{\partial q_D}\right) = E.\] \(E\) 为哈-雅方程解得常数 \(\alpha_1\), 表明 \(H\) 为运动常数.

3.2.3. 可分离系统

  • 若体系得哈雅方程可以用分离变量求解, 则称之为可分离体系.
  • 若所有变量均可分离,则称完全可分离体系.

是否可分离还与广义坐标的选取有关.

对于不含时哈雅方程: \[H \left( q, \frac{\partial W}{\partial q} \right) = E.\] 若满足下列特殊形式: 方程左式某一坐标 \(q_l\) 和与之相对应偏导数 \(\partial_{q_l}W\), 总是以某种函数形式 \(\phi(q_l, \partial_{q_l}W)\) 出现且 \(\phi_1\) 不以其他变量做自变量.

则不含时哈-雅方程可以分离出 \(q_l\) 变量的函数, 即 \[H \left[ q_{[l]}, \frac{\partial W}{\partial q_{[l]}}, \phi \left( q_l, \frac{\partial W}{\partial q_l} \right) \right] - E = 0.\] 可得 \[W = W_{[l]}(q^{[l]}) + W_l(q_l).\] 此时,哈-雅方程化为 \[H \left[ q_{[l]}, \frac{\partial W_{[l]}}{\partial q_{[l]}}, \phi \left( q_l, \frac{\mathrm{d} W_l}{\mathrm{d}q_l} \right)\right] - E = 0.\] 可解得 \[\phi_1 \left( q_l, \frac{\mathrm{d}W_l}{\mathrm{d}q_l} \right) = f \left( q_{[l]}, \frac{\partial W_{[l]}}{\partial q_{[l]}} \right).\] 与 \(q_l, \mathrm{d}_{q_l}W_l\) 无关. 可得: \[\phi_1 \left( q_l \frac{\mathrm{d} W_l}{\mathrm{d} q_l} \right) = C_2 = \text{const}, \;H \left[ q_{[l]}, \frac{\partial W^{[l]}}{\partial q_{[l]}}, C_2 \right] = E.\] 对于完全可分离体系,哈雅方程可分离出\(D\)个常微分方程.

当某个坐标为循环坐标时, \(H\) 不显含该坐标, \(\phi_1\) 退化为 \[\phi \left( \frac{\mathrm{d}W_l}{\mathrm{d}q_l} \right) = C_2 \Rightarrow \frac{\mathrm{d}W_l}{\mathrm{d}q_l} = C_2^{\prime} \Rightarrow W_l = C_2^{\prime}q_l.\]

将求解哈密顿方程变换成求解哈-雅方程不一定使求解更容易.

4. 作用量变量与角变量

4.1. 多周期体系与作用量变量

完全可分离体系

一个哈密顿量不显含时间的体系, 若存在一组广义坐标 \(q_1,\dots, q_D\), 使得其不含时哈-雅方程 \[H \left( q, \frac{\partial W}{\partial q} \right) = E.\] 存在以下形式解: \[W = W(q_1, q_2, \dots, q_D, \alpha_1, \alpha_2, \dots, \alpha_D) = W_k(q^k, \alpha_1, \dots, \alpha_D).\] 则称该体系为 完全可分离体系.

对于完全可分离体系, \(p_k = \partial_{q_k}W_k = p_k(q_k, \alpha)\).

多周期体系

完全可分离 哈密顿体系投影在每个 \(q_k, p_k\) 相平面,

  1. \(p_k = p_k(q_k, \alpha)\) 闭合, 此事称为 天平动.
  2. \(p_k(q_k, \alpha) = p_k \left[ \left( q_k + q_k^T \right), \alpha \right]\), 但 \(q_k\) 不一定为时间的周期函数, 此时称为 转动.

对所有\(k\)都满足以上两性质,则称之为 多周期体系.

4.1.1. 多周期体系利用周期性简化问题

\[J_k = \frac{1}{2\pi}\oint p_k \mathrm{d}q_k,\;k=1, 2, \dots, D.\] 以上一组量为常量,称为作用量. (其量纲为作用量量纲) 考虑到 \(p_k = p_k(q_k, \alpha)\), \(J_k\) 仅与 \(\alpha\) 有关, 与 \(q_k\) 无关. \[J_k = J_k(\alpha)\Rightarrow\alpha_k = \alpha_k(J).\] 哈密顿特征函数有: \[W(q, \alpha) = W_k(q^k, \alpha_1,\dots,\alpha_D) = \tilde{W}_k(q_k, J_1,\dots,J_D).\] 认定 \(J\) 为变换后广义动量 \(P_k\). 取哈密顿特征函数为生成函数, 生成变换仍是特征变换. 经其变换, \(K = H + \partial_t W = H = E = \text{const}\). 记正则动量 \(P = J\) 的共轭坐标 \(Q = \phi\), 由生成函数可得: \[\phi_k = \frac{\partial \tilde{W}}{\partial P_k} = \frac{\partial \tilde{W}(q, J)}{\partial J_k}.\] 称为 角变量. 易证\(\phi\)是无量纲量. 另外有: \[\dot{\phi}_k = \omega_k = \text{const}, \; \dot{P}_k = 0.\] \(\omega_k\) 对时间为常数, 依赖于 \(J\): \(\omega_k = \omega_k(J)\).

4.2. 运动周期

角变量\(\phi_k\) 在一个周期的变化 \[\Delta\phi_k = \oint \mathrm{d}\phi_k = 2\pi.\] 设 \(\Delta \phi_k = \omega_k\tau_k\), \(\tau_k\) 为 \(q_k\) 运动一周期时间. 可得: \[\omega_k = \phi_k = \frac{\partial K}{\partial P_k} = \frac{\partial E(J_1, \dots, J_D)}{\partial J_k}.\]

4.3. 简并

不同相平面周期不尽相同.

简单周期体系
各相平面周期之比为有理数.
条件周期体系
有两个或更多相平面内运动周期之比不是有理数.

条件周期体系运动状态永远不能完全复原,简单周期函数反之.

简单周期体系有: \(n_1\omega_1 = n_2\omega_2 = \cdots = n_{D-1}\omega_{D-1} = n_D\omega_D\). 因此存在 \(D-1\) 个以下形式独立关系式: \[n_k^{(l)}\omega^k = 0,\;n_k^{(l)}\in \mathbb{N},\;l=1, \dots, D-1.\]

简单周期体系完全简并. 简并条件可用于减少作用量数目.

\(m\)重简并体系
如果一个体系存在 \(m\) 个以上形式的关系时, 称其为 \(m\) 重简并体系.

4.3.1. 利用简并减少系统作用量数目

做正则变换 \((\phi, J)\rightarrow (\tilde{\phi}, \tilde{J})\): 其生成函数为: \[U_2(\phi, \tilde{J}) = \tilde{J}_l n_j^{(l)}\phi^j + \tilde{J}_l\phi^l\]. 第二项实际上对 \((\phi_k, J_k)\), \(k = m+1, \dots, D\) 做恒等变换. 得到:

\begin{equation*} \tilde{\phi}_l = \begin{cases} n_j^{(l)}\phi_j&, l = 1,\dots, m\\ \phi_l &, l = m + 1, \dots, D \end{cases} \Rightarrow \tilde{\omega} = \begin{cases} 0 &, l = 1, \dots, m\\ \omega_l &, l = m+1, \dots, D \end{cases} \end{equation*}

即对于\(m\)重简并体系, 经过正则变换, 只剩下 \(D-m\) 馆不同的频率, 自由度降低.

日期: 2019-11-10 Sun 11:01

作者: yuandi